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ABSTRACT 

We calculate the Yukawa corrections of order 0(a ew m^, b \/rn^ v ) to charged Higgs 
boson production in association with a top quark at the Tevatron and the LHC The 
corrections are not very sensitive to the mass of the charged Higgs boson and can 
exceed —20% for low values of tan j3, where the contribution of the top quark is large, 
and high values of tan/3 where the contribution of the bottom quark becomes large. 
These Yukawa corrections could be significant for charged Higgs boson discovery 
searches based on this production process, particularly at the LHC where the cross 
section is relatively large 



PACS number: 14.80.Bn, 14.80. Cp, 13.85. QK, 12.60. Jv 



1. Introduction 



There has been a great deal of interest in the charged Higgs bosons appearing in 
the two-Higgs-doublet models(THDM)[l], particularly the minimal supersymmetric 
standard model(MSSM)[2], which predicts the existence of three neutral and two 
charge Higgs bosons ho,H,A, and H ± . The lightest neutral Higgs boson may be 
difficult to distinguish from the neutral Higgs boson of the standard model(SM), but 
charged Higgs bosons carry a distinctive signature of the Higgs sector in the THDM 
and MSSM. Therefore, the search for charged Higgs bosons is very important for 
probing the Higgs sector of the THDM and MSSM and, therefore, will be one of 
the prime objectives of the CERN Large Hadron Collider (LHC). At the LHC the 
integrated luminosity is expected to reach L = per year. Recently, several 

studies of charged Higgs boson production at hadron colliders have appeared in the 
literature[3,4,5]. For a relatively light charged Higgs boson, m H ± < m t — m b , the 
dominate production processes at the LHC are gg — > tt and qq — > ti followed by the 
decay sequence t — > bH^ 1 — ► 6r ± z/ r [6]. For a heavier charged Higgs boson the dominate 
production process is gb — > tH~ [7,8,9]. Previous studies showed that the search for 
heavy charged Higgs bosons with m#+ > m t + m b at a hadron collider is seriously 
complicated by QCD backgrounds due to processes such as gb — > tib, gb — > tib, and 
gg — > tibb, as well as others process [8]. However, recent analyses [10, 11] indicate 
that the decay mode H + — > r + u provides an excellent signature for a heavy charged 
Higgs boson in searches at the LHC. The discovery region for H ± is far greater than 
had been thought for a large range of the (m#±, tan/3) parameter space, extending 
beyond m#± ~ ITeV and down to at least tan/3 ~ 3, and potentially to tan/3 ~ 1.5, 
assuming the latest results for the SM parameters and parton distribution functions 
as well as using kinematic selection techniques and the tau polarization analysis [11]. 

The one-loop radiative corrections to H~t associated production have not been 
calculated, although this production process has been studied extensively at tree- 
level[7,8,9]. In this paper we present the calculations of the Yukawa corrections to 
this associated H~t production process at both the Fermilab Tevatron and the LHC 



1 



in the THDM. These corrections arise from the virtual effects of the third family(top 
and bottom) quarks, the charged and neutral Higgs bosons, as well as the Goldstone 
bosons. The one-loop QCD corrections are probably more important, but are also 
more difficult to calculate, and we will present these calculations in a future publica- 
tion^]. 



2. Calculations 

The tree-level amplitude for gb — > tH~ is 

M = Mi s) + M^\ (1) 

where and M^p represent the amplitudes arising from diagrams in Fig. 1(a) and 
Fig. 1(6), respectively. Explicitly, 

M (s) = W9s -u{p t ) [2m t cot (3$P L + 2m b tan (3p^P R - m t cot (3^ fiP L 

V2m w {s — m() 

-m b ttmPY /^PhMp^e^T", (2) 

and 

= W9s -u(p t ) [2m t cot (3pt P L + 2m b tan (3p?P R - m t cot 0-f fcP L 

V 2m w {t — mf) 

-m b twP>y'' / kP R ]u(p b )e lt (k)T? j , (3) 

where T a are the SU (3) color matrices and s and i are the subprocess Mandelstam 
variables defined by 

s = (p b + k) 2 = (pt+PH-) 2 , 

and 

t = (p t - kf = (p H - - p b f. 

Here the Cabbibo-Kobayashi-Maskawa matrix element VcKM[bt] has been taken to 
be unity. 
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The Yukawa corrections of order 0{a ew m 2 , b Jm 2 v ) to the process gb — > H~t arise 
from the Feynman diagrams shown in Figs.l(c)-l(v) and Fig. 2. We carried out the 
calculation in the t'Hooft-Feynman gauge and used dimensional regularization to 
regulate all the ultraviolet divergences in the virtual loop corrections using the on- 
mass-shell renormalization scheme[13], in which the fine-structure constant a ew and 
physical masses are chosen to be the renormalized parameters, and finite parts of the 
conterterms are fixed by the renormalization conditions. The coupling constant g is 
related to the input parameters e, m w , and m z by g 2 = e 2 /s 2 u and = 1 — m 2 u /m 2 z . 
The paramerter j3 in the THDM we are considering must also be renormalized. Follow- 
ing the analysis of ref . [14] , this renormalization constant was fixed by the requirement 
that the on-mass-shell H + lui coupling remains of the same form as in Eq.(2) of ref. [14] 
to all orders of perturbation theory. Taking into account the 0{a ew m 2 , b Jm 2 v ) Yukawa 
corrections, the renormalized amplitude for the process gb — > tH~ can be written as 



where 8M V ^ S \ 8M V ^\ 8M< S \ 8M"®,8M V *W, 6M V *®, 8M b( ~ s \ and 8M b ^ represent 
the corrections to the tree diagrams arising, respectively, from the gbb vertex di- 
agram Fig. 1(c), the gtt vertex diagram Fig. 1(e), the bottom quark self-energy di- 
agram Fig. 1(g), the top quark self-energy diagram Fig.l(i), the btH~ vertex dia- 
grams Figs.l(k)-l(m) and Figs.l(o)-l(q), including their corresponding counterterms 
Fig. 1(d), Fig. 1(f), Fig. 1(h), Fig.l(j), Fig.l(n), and Fig.l(r), and the box diagrams 
Figs.l(s) — 1(f). J2i 3M l then represents the sum of the contributions to the Yukawa 
correctons from all the diagrams in Figs.l(c)-l(v). The explicit form of 8M l can be 
expressed as 



where the C l are coefficients that depend on s, t, and the masses, and the f\ are 



ren 



M (s) + M W + 8M Vl(s) + 8M Vl(t) + 8M s{s) + 8M s(t) + 5M V2(s) 
+5M V2{t) + 8M b{s) + 8M m = M,S s) + M W + £ 8M l , 



(4) 



8M 1 
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form factors; both the coefficients C l and the form factors f\ are given explicitly in 
Appendix A. The corresponding amplitude squared is 

£|M ren | 2 = £|M (s) + M^ 2 + 2i?e£[(£ 5M l )(M^ + M«)t], (6) 

z 

where 

£|M^ + M W | 2 = jrp^-i f. 1 2 . 2 [(m 2 cot 2 /3 + m 2 tan 2 /3)(p 6 -A;p t -fe-m 2 p t -A; 

,2„ „ \ , r. 2 2/ . 2\ 



+ 2p b • A;p 6 • pt - m fe p 6 • p t ) + 2m b m t (p b ■ k - m b )] 
1 

+ 2m 2 b m 2 t {p t ■ k - m 2 )] + 



+ — — —J^i( m t cot2 P + m t tan2 P)(Pb • kp t • k + m 2 p b ■ k - m 2 p b ■ p t ) 

1 



(s-m 2 )(t-m 2 ) 
x [(m 2 cot 2 /3 + m 2 b tan 2 /3) (2p fe • kp t ■ k + 2p b ■ kp b ■ p t - 2(p b ■ p t ) 2 

- m 2 b p t - k + m 2 p b - k) + 2m 2 b m 2 (p t - k-p b - k-2p b -p t )]}, (7) 



A„ 12 



and 



,4„ 12 



5jSM*(MW)t = - //% ^^/if/' (9) 

Here the color factor N c = 3 and /if ^ and /if ^ are scalar functions whose explicit 
expressions are given in Appendix B. 

The cross section for the process gb — > tH~ is 

1 



-max , 



«r = / ^\M ren \ 2 dt (10) 



with 

^2 j ™2 



*min = ^ +m 2 H - \/(s- (m t + m H -) 2 )(s - (m t - m H -) 2 )/2, 

and 



_|_ ffl^ 1 — Q j 

tmax = — + J(s-(m t + m H -) 2 )(s-{m t -m H -) 2 )/2. 
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The total hadronic cross section for pp — > gb — > ii? can be obtained by folding the 
subprocess cross section a with the parton luminosity: 

a(s) = f 1 dz^a(gb^tH~ at s = z 2 s). (11) 

Here ^/s and are the CM energies of the pp and #6 states , respectively, and dL/dz 
is the parton luminosity, defined as 

Jz =2Z L ^fv/P&rifg/P^/B't*)' ( 12 ) 
where f q /p(x,/j,) and f g /p(z 2 /x,fj,) are the quark and gluon parton distribution func- 
tions. 

3. Numerical results and conclusion 

In the following we present some numerical results for charged Higgs boson pro- 
duction in association with a top quark at both the Tevatron and the LHC. In 
our numerical calculations the SM parameters were taken to be m w = 80.33C7e^, 
m z = 91.187GeV, m t = 176GeV, m b = 4.5GeV , a s = 0.118, and a ew = ±. For sim- 
plicity, we also used the relations from the MSSM between the Higgs boson masses 
m h ,H,A,H ± an d the parameters a and f3, and chose m^± and tan/? as the two in- 
dependent input parameters. And we used the CTEQ5M[15] parton distributions 
throughout the calculations. 

Figures 3(a) and 4(a) show the tree-level total cross sections as a function of 
the charged Higgs boson mass for three representative values of tan (3. For m H ± = 
200GeV the total cross sections at the Tevatron are at most only a few fb for tan (3 = 
2, 10, and 30, and for m H ± = 300GeV the total cross sections are smaller than 1 
fb for all three values of tan/3. However, at the LHC the total cross sections are 
much larger: the order of thousands of fb for m H ± in the range 100 to 300GeV and 
tan (3 = 2 and 30; and they are hundreds of fb for the intermediate value tan (3 = 10. 
For low tan f3 the top quark contribution is enhanced while for high tan f3 the bottom 
quark contribution becomes large. These results agree with ref.[8,9] and, it should be 
noted, are larger than the W ± H ± associated production cross section at the LHC [4]. 
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In Figs. 3(b) and 4(b) we show the corrections to the total cross sections relative to 
the tree-level values as a function of m H ± for tan (3 = 2, 10, and 30. These corrections 
decrease the total cross sections significantly for a wide range of the charged Higgs 
boson mass, especially for the smaller values of tan (3 where the top quark contribution 
is greatly enhanced. In particular, for tan/3 = 2 the corrections exceed —20% for 
m H ± below 300Ge\/ and reach more than —25% for m H ± below 200GeV at both the 
Tevatron and the LHC. 

In conclusion, we have calculated the Yukawa corrections of order 0(a ew m^/m^ v ) 
to the cross section for charged Higgs boson production in association with a top quark 
at the Tevatron and the LHC. These corrections decrease the cross section and are 
not very sensitive to the mass of the charged Higgs boson, but depend more strongly 
on tan (3. At low tan/3 the top quark contribution is enhanced while at high tan/3 
the bottom quark contribution becomes large. For m#± in the range 100 to 300 GeV 
the Yukawa corrections are as large as —30% for tan (3 — 2, then become smaller for 
the intermediate value tan (3 = 10, but increase to nearly —20% for tan/3 = 30. 

This work was supported in part by the National Natural Science Foundation of 
China, the Doctoral Program Foundation of Higher Education of China, the Post 
Doctoral Foundation of China, a grant from the State Commission of Science and 
Technology of China, and the U.S. Department of Energy, Division of High Energy 
Physics, under Grant No.DE-FG02-91-ER4086. S.H. Zhu also gratefully acknowledges 
the support of the K.C. Wong Education Foundation of Hong Kong. 
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Appendix A 



The coefficients C l and form factors f\ are the following: 

2 2 2 

(jVi(s) _ m b (jVxit) _ m t qs(s) _ "'i- 



m 2 w (s-m 2 b y m 2 w {i-m 2 t y m^(s-m^) 2 ' 

(js{t) _ m t (jV 2 (s) _ 1 (jV 2 (t) _ ' 



fjb(s) = Qb(t) 



m 2 v {i-m 2 ) 2 '' s-m 2 ' i - m 2 ' 

1 



m w 



V W [m b (92 l(s) -93 l(s) )-2p b -kx 9 ^% 

V {2) [rn b (9^ {s) -9^ s) )-2p b -kxg^% 

„<W lW + 9l l{s) ) + m b ( 9 ? W + £ {S) ) + 2 Pb ■ k x gW], 

U ( W lW + g? W ) + m b (gl l(s) + <fc* W ) + 2 Pb ■ k x g™], 

^[-(/ 1(S) +^ (S) )+- 6 (^ (S) + ^ (S) )], 
V { %: i{s) + 2g^ + m b (g^-g^)], 

V { %^ s) + 2g^ s) +M9^ s) -9^ s) )], 
2 Pb -kg^\ f^ = 2p b .kg^\ 

2g\ 2is) + 2m t cot /3SA^ H ~ - 2m t g^ (s) + 2m b gf {s \ 
2g^ s) + 2m 6 tan/MA*f " + 2m b gf {s) - 2m t g\^\ 

_1/V 2 (s) f V 2 (s) _ _1 f V 2 (s) 

2^3 i J8 ~ ' 

2 Pt -kg^\ f?® = 2p t .k g ™ 

2g™ + 2m, cot /36A^ H ~ - 2m t gf {t) + 2m b gf {t) , 

2gl 2{t) + 2m 6 tan/WAg ff ~ + 2m b g? {t) - 2m t g V ^\ 

1 ,v 2 (t) ,v 2 (t) _1 ,v 2 (t) 

2 J 5 ) J 8 2 •'6 5 

2^ 2 )p b .A;[yi (s) + m 6 (^ ) + ^ s) )], 
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/ 3 S(S) = 2^[2m b gf s) + 2{m\ + Pb ■ k)g s 2 ^ + (m 2 b + 2 Pb • k)g s 3 {s) + m 2 gf\ 
/ 4 S(S) = 2^[2m b gt {s) + 2{m\ + p b • k)g s 2 {s) + m 2 b g^ s) + (m 2 b + 2p b • k)gt% 

rs(s) _ 1 ns(s) f s(s) _ 1 »s(s) 

+2p t • k(D 13 + D 26 ) + 2p b -p t {D 13 + D 2 5)](-Pb, -k^p^m^m^m^mj) 
+— E d 3) {^ (2) [2m 6 (-3L> 3 i2 + (1 - 0)^27) + m?(A, + D 12 

m W i=H°,h°,G°,A° 

-D 22 - D 32 ) - m 2 m b (D 23 + 2D 39 ) - 2m b p b ■ k(2D 36 + D 2A + d(D + D 12 )) 

+2m bPt ■ k(D 25 + D 310 ) + 2m bPb ■ p t (D 26 + 2D 38 )} + r]^[2m t (-3D 313 

+ (1 + QD 27 ) - m 3 t (D 33 + (1 + (i)D 23 ) + m 2 b m t (D 13 - 2D 38 + (1 + Q(D 

-Dm)) + 2m tPb ■ k(D 13 - D 3W - (1 + Q)(D 12 + D 24 )) + 2m tPt ■ k(2D 37 

+ (1 + Ci)D2s) + 2m t p b ■ p t (2D 39 + (1 + d)D 26 )]}(-k, -p b ,p t , m b , m b , m h m t ), 

fs (S) = E §dMd) 2m ^ D " + D 21 + (1 + Ci)(D + D n )} - V^m^ + D 25 )} 

(id) 

(-p b ,-k,p t ,m i ,m b ,'m b ,'m j ) + H^± ^ {r] (1) [-AD 27 + 2m 2 b (D 22 

m W i=H ° ,h° ,G° ,A° 

-D - (1 - Q)(D 12 + D 22 )) + 2m 2 (D 23 - (1 + Q)D 26 ) + 4p t • k(D 26 
-D 25 )] +ri {2) 2m t m b (l + Q( D 22 ~ D 12 - D 26 )}(-k, -p b ,p t ,m b ,m b ,mi,m t ), 

/5 (S) = E £(S)^m')(- 2 ^) P» + (1 + CO As] + ^ ( ( 3)2m tJ D 23 } 

(-p b , -k, Pt , rrn, m b , m b , m 3 ) + ^ E d 3) {^ (1) [12 Aos 

m W i =H t h0,G ,A° 

+2m 2 b (2D 38 - D 13 + (1 - Ci)(£>i 3 + D 26 )) + 2m 2 (D 33 + (1 + Q)D 23 ) 
+4p b ■ k(D 25 + D 3W ) - Ap t ■ k{D 23 + 2D 37 ) - Ap t ■ p b (D 23 + 2D 39 )] 
+ri (2) 2m t m b (l + d)(D 13 + £> 2 3 - D 26 )}(-k, -p b ,pt,m b ,m b ,mi,m t ), 

f Hs) = # VW^g) ^3))> 

fr (S) = E^Ml)(-rn b )[D 11 + (l + Ci)D } + 

(id) 
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(-p b ,-k,pt,mi,m b ,m b ,mj)-\ — E d 3) {^ (1) [6(^27 

m W i=H °,h ,G ,A° 

-D 311 ) + m 2 b (D u - 2D l2 - 2D 22 - 2D m + (1 + Q{D Q + D 12 )) - m\{2D 23 
+2D 31 + (1 + d)D 13 ) - 2p b ■ k(D 12 + 2D 24 + 2D U ) + 2p t ■ k(D 13 + 2D 25 
+2D 35 ) + 2p t ■ p b (D 13 + 2D 26 + D 3W )} + r]^m t m b (l + Q)(D 12 - D 13 
-D )}(-k, -p b ,p t ,m b ,m b ,mi,mt), 

f b s (s) = f? S \v m ^rj {2 \vZ^& 

/9 (S) = E^i)^) 2 ^ 12 + £) 24)](-P6,-A:,Pt, 77^, 7716, 77l 6 , mj) 

+ — E d 3) {^ (1) 2m t [-£>i3 - ^26 + (1 + 0)(D 12 + L> 24 )] 

m ^ i=H°,h ,G°,A° 

+7] {2) 2m b [-D 22 + L> 2 4 + C*(A> + 2 ^i2 + D 24 ]}(-k, -p b ,p t ,m b ,m b ,mi,m t ), 

/5 s) = / 9 6 % (1 w^iWS))> 

/n S) = E ^(U)' - ^') 2 ^ 13 + ^26)](-P6, -k,p t , mi, m b , m b , rrij) 
+ rn f n b E e (3) {77 (D 2mt[jD23 _ (1 + Ci)jD25] 

m W i=HO,hO,G°,A° 

-r]^ 2) 2m b [-D 26 + D 25 + d(D 13 + D 25 \}(-k, -p b ,p t ,m b ,m b ,mi,m t ), 



Ms) 
J 12 


= A?^- 


^ (2) ^(S)^^ ( 3))' 










f m) 


= f? l{s \u), 




f Vi(t) _ 
J5 


fZ l(s \u), 


fVi(t) 






= fs l(s) (U), 




fVi(t) _ 
ill _ 


-f™(U), 


fVi(t) _ 
J 12 






= A s(s V), 




Mt) _ 
J5 — 


ff s \u), 


Mt) _ 
J6 — 


fs (s \u), 


f S { t) 


= -\ft\u), 




jb{t) 


= f b(s \u), 


Mt) _ 

J2 — 


f b 2 (S \U), 


Mt) 
J3 


= tf'\U), 


f b{t) = f b{s \u), 


Mt) 

J5 — 


f b(s \u), 


Mt) 
J6 — 


f b 3 (S \U), 


f b { t) 


= fs is) (U), 


fs {t) = f b 7 is \u), 


Mt) _ 
J9 — 


-AtV), 


Mt) _ 
J 10 — 


-AtV), 


Mt) 
ill 


= -f b 9 (S \U), 


fi? = -fi { o\u), 











Here the sums over run over (H°,H~), (h°,H-), (H°,G~), (h°,G~)and(A°,G-) 
and U is a transformation defined by 
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m t ^m b , 77^^77(2), vfy^vffj), £(3)^(3)' 

and D , Dij, Dijk are the four-point Feynman integrals [16]. All other form factors 
fl not listed above vanish. In the above expressions we have used the following 
definitions: 



,( 2 ) 



t(2) 
t(3) 

t(3) 
Sff- 



_„(!) 

— 'l(h°,H-) 


= 77« = 


777,5 tan /3, 


_„(2) 


= t/ 2 ) = 


m t cot /5", 


cos 2 a 




sin 2 a 


cos 2 p 1 ' 


cos 2 f3 ' 


sin 2 a 


t(2) 


cos 2 ct 


sin 2 p 1 ' 


sin 2 (5 ' 



t(3) 



m 



sin a cos a 



sin /3 cos /3 ' 



— | cot^ 



mt 



= tan 2 /?, 



4 S » = 1. 



(i) _ (i) 



„(2) 

'l(H°,G- 

£$ = tan 2 /5\ 
cot 2 /3, 



„(2) 



t(2) 



t(3) 



t(3) 



(2) -^|tan 2 /5, 



Si? 

d 4 ^ = cot 2 A 



(1) 

^(AO,G-) — m& > 
(2) 



<r(i) 
^G° 



t(2) 
SG° 



1, 
1, 



t(2) _ 
^G - 2 ' 



$(H°,H- 

Ai) 

$(h°,H- 

Hh°,G~ 

?(AO,G- 
t(2) 

$(H°,H- 
t(2) 

t(2) 

?(if°,G- 

t(2) 
?(hf,G- 

t(2) 

^(A°,G- 



cosa r . _ . mz ^ . ^ ., 
2m fe mvj/ cosfp — a) — cos 2p cosfp + a) , 

cosp 2 COSfe'iy 

sin a . . 

—2mb t: m w smfp — aH — cos 2p smfp + a)\ , 

cos p 2 cos 6V 

cos a , rn - 



'~ r~ vv 

m b tSj^w sin(/3 — a) f— sin 2/3 cos(/3 + a)] 

cosp cosfcV 



sin a 



[mvK cos(/3 — a) ^— sin 2/3 sin(/3 + a)] , 

COS 6% 



cos p 1 
m^ra^ tan p 1 , 

sina r . . mz n „ , n xn 
2m t 7:\mw cosfp — a) — cos2p cosfp + a) \, 

sin p 2 cos tV 

cosa; r . . . . . ., 

2m t 7:[mw smfp — aH — cos 2p smfp + a) , 

ship 2costV 

sin a m z . 

ra* 1 smfp — a) — sin 2p cosfp + a J , 

snip cosfcV 
cosa r . . mz . nn . , n Nn 

— m t 7:[mw cosfp — a) — sm2p smfp + ml, 

sin p cos 9w 

m t mw cot P, 
Cm = (h° = (h- = ~(ao = -Cg" = ~Cg- = 1, 
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E tPiil - 2C 24 + m 2 b (-2Cu + Ci 2 - C 21 + C 23 ) - s(C 12 + C 23 )] 

-k,mi,m b ,m b ) + [-F + F 1 + 2m 2 b G 1 - (1 + C^m^cKm 2 , m 8 , m b )}, 
E 2{^ (1) [i - 2C 24 + m 2 C + m 2 (-G) - 2C n + C 12 - C 21 + C 23 ) 

i=H-,G~ 

-s{C 12 + C 23 )](-p b , -^m l) m t) m t ) + [^(-Fq + F x ) - 2m t 2 CiGo 
+ C (3) )(Gi - 0G )](m 2 ,m 4 ,m t )}, 

^ (s) (#-d 3) ;^^-,G-), 

E & ] 2m b [C + 2C n + C 21 + Ci(C + C n )](-p 6 , -fc,m m b ,m b ) 

i=H ,h°,G°,A 

.( 3 ),„ _ „ s m 2 



+ ]T W£i (Co + 2C n + C 21 ) + -4Ci(C + C n )](-p 6 , —k : m i: m t: m t ) : 

i=H-,G~ m b 



E d 1)m &(Co + Cn + OC'o)(-p fe , -fc, mi, m b ,m b ) 

i=H°,h°,G a ,A 

- E 2m fe [^ 3) (C + C n ) H |C*Co](-Pb, -k,mi,m t ,m t ), 

i=H-,G~ m b 

E ^ 1} 2(Ci 2 + C 23 )(-p 6 , -k,mi,m b ,m b ) 

i=H ,h°,G°,A 

+ E ^f^C^ + C 24 )(-P6, -k,mi,m t ,m t ), 

i=H-,G~ 

E ^d 3) {^ {1) [-4 + 4^24 + m 2 (C + 2C n + 0(C + C u ) + C 21 

i=ff°,/i°,G°,A0 m ^ ^ 



-C12 - C 23 ) + (C 22 - C 23 ) + s(C 12 + C 23 )] + ?? (2) m 6 WCo + Ci(C 
+Cn)]}(-p t , -p H -,mi, m t , m b ) + E {^(3) ? ?(S) m *[( 1 + CO Co + C i2 ] 

(-PH-,~Pt, m^rrii, m t ) + £(^)fa(?j) m t(Co + C i2 ) + ^m^Co] 
(-PH-:-Pt:m,m j: m b )}, 

E ^d 3) {^ (1) ^[Co + C n + C*(Cb + C 12 )] + f/ 2 >C^C 12 } 
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(~Pt, -PH-,mi, m t , m b ) + E — [Avu]) (d + C n ) 

{ij) m w ^ ^ 

(~Ph- , ~Pt, mj,mi, m t ) + f [v$j)(C + Cn) , -p t , m, rrij, m b )\, 

E ^d 3) {^ (1) [-^ + 4C 24 + m 2 (q, + 2C n + Q(C + C n ) 

+C21 — C12 — C23) + J^h- (C22 — C23) + t(C\2 + C23)] 
+r] ( - 2) m b m t [C + 0(C + C u )]}(-p b ,p H - ,mi,m b ,m t ) 

+ E ^{^( ( 3)^(S) m *[( 1 + Ci)Co + Ci 2 ](-Ph- ,p b ,mj,mi,m b ) 
mw 

+£(5)fo(ii m6 ( C ° + Cn) + ri$j) m t(iC }(-PH-,Pb, m, rrij, mt)}, 

E -^d 3) {?7 (1 VO, + C n + 0(C + C 12 )] + ^Cim*C 12 } 

i=H ,h°,G°,A m W 

(-p b ,PH-,mi,m b ,m t ) - E —[t{i])V(i))( C o + Cii)(-pjr-,p 6 , rrij, m*, m t ) 

mw 

+£(5)^(S)( C ° + Cn)(-p ff -,p 6 ,m i ,m : ,-,m t )], 

3 Vi % (1 w% ( 3) 

E ^Cf^-C^oU + k, mi, m b ) + [C,F - 2m 2 (l + 0)C 

i=HO,h° ,G<-\A° 

2 

+2m^G 1 ](m^m ! ,m () )}+ E 2m &{ |0^o(p& + fc, ™i, "it) 

i=H-,G~ m *> 

2 

+ [-2m, 2 0G + m2(e i (1) +d 3) )(Cx -OCo) + cA i71 o]K,m J ,m,)}, 
E d'V^o + ^O^ + ^m^m,), 

E dVo - F x - 2mgGi + 2(1 + OKC ]K,™^) 

+ E 2{C (1) (-^o + i 71 i)(P6 + A;,m l) m t )-[d 1) (-^o + i 71 i) 

i=H-,G~ 

-20m, 2 G + m^(d 1) +C (3) )(Ci - C,C )](m 2 , m*, m t )}, 
^-(1 - cot 2 9 w )[2m 2 t (ln ^ - 1) + m 2 + m 2 - + m 2 F 
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+(m 2 b -m\- 2m 2 w )F 1 ](m 2 w ,m b ,m t ) + cot 2 6 w {-\[{g b v ) 2 + (g b A ) 2 



Hg'v? + (<&)>! + Wv? + (^) 2 )(2m 2 In ^ + m 2 F - 2m 2 z F 1 ) 

"((^) 2 - (^) 2 )3m 2 F ](m|,m t ,m i ) + [((<&) 2 + (^) 2 )(2m 2 In ^ + 

/i 

-2m|F0 - ((<&) 2 - (^) 2 )3m 2 F ](m|,m 6 ,m 6 )} + ^[(cot 2 /? - l)m 2 F 
+ (m 2 -ra\- 2m 2 cot 2 + (m 2 cot 2 /3 + m\ tan 2 /3 + 2m 2 b )m 2 G 
— (m 2 cot 2 /3 + m 2 tan 2 P)m 2 H -Gi](m 2 H - , m t , m&) 

+ E A-W&Vi - - 2m 2 (l + Ci)G + 2m 2 G 1 ](m 2 , m„ m 6 ) 

i=H°,h°,G°,A° lm W 

-m 2 t g\-F Q + F 1 - 2QF + 2m 2 (l + Q)G - 2m 2 G 1 ](m 2 ,m l ,m t )} 
+ E ^W 2 [d 1) (-^o + Fi)-2m 2 OGo + m 2 (d 1) +d 3) ) 

i=H-,G- m W 

x(Gi - CiGoKm^m^mt) - m 2 [-^f C^o + d 2) (-^o + F) + 2m 2 0G 

"it 

-m 2 (d 2) +d 4) )(Ci -CiGoMmlm^m,)}, 
^r H ~ = E ^{m 2 d 2) [-F + F 1 -2m 2 (l + 0)G + 2m 2 G 1 ](m 2 ,m l ,m t ) 

i=H°,h°,G°,A° Am W 

-m 2 C (1) [-F + Fx - 2QF + 2m 2 (l + QG - 2m 2 G 1 ](m 2 , m,, m fe )} 

+ E ^K[d 2 V^o + Fi)-2m 2 OG + m 2 (d 2) +d 4) (Gi-OG ] 

i=H-,G~ m W 

{ml m., m b ) - m \\-^_QF, + £\-F + F,) + 2m 2 QG - m 2 b (^ + &< s >) 
x(Gi - CiGo)](m 2 b ,mi,m t )}. 

Here C , are the three-point Feynman integrals [16] and C 24 = — + C 2 4, while 

.-g 2 y(l - y) + m 2 (l - y) + m|y 



A* 2 



F n (q,m 1 ,m 2 ) = f dyy n \n[- 

G„( 9 ,m„m 2 ) = -^*^ wl _„ )+nh » (1 _ B)+n4) 



2„ ,' 



and 



^y = ^-^ sin2 ^, #a = ^, ^ = -^ + ^sin 2 ^, 9a = ~\, 
which are the SM couplings of the top and bottom quarks to the Z boson. 
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Appendix B 



hi = Ami cot (3(2p b ■ k — ml) — Ami tan/3(p& • p t + pt • k), 

h^ = — Am b m t cot (3 (p b ■ p t + p t ■ k) + Am b m t ta,n (3 (2p b ■ k — ml), 

h^ = 2m t cot (3(2p b ■ kp b ■ p t - m 2 b p t ■ k - 2m 2 b p b ■ p t ) + 2m 2 b m t tan/3(p 6 • k - 2ml), 

= 2m\m b cot f3(p b ■ k - 2m 2 ) + 2m b tan/3(2p 6 • kp b ■ p t - m\p t ■ k - 2m 2 b p b ■ p t ), 

h\° ] = 2m t cot(3(m 2 p b ■ k -2(p b ■ p t ) 2 ) + 2mlm t tan(3(p t ■ k -2p b ■ p t ), 

h { 6 s) = 2m 2 t m b cot (3(p t ■ k - 2p b ■ p t ) + 2m b tnn (3(m 2 t p b ■ k - 2(p b -p t ) 2 ), 

= Am t cot f3{mlp t ■ k — 2p b ■ kp b ■ p t — 2p b ■ kp t ■ k) — Am\m t tan j3p b ■ k 

h^ = —Am 2 t m b cot (3p b ■ k + Am b tan (3(mlp t • k — 2p b ■ kp b ■ p t — 2p b ■ kp t ■ k), 

hg = Ami c °t ftPb ' k(p b ■ k — ml) — Ami tan (3p t ■ k, 

hy] = — Am\m t cot (3p t ■ k + Am b m t tan (3p b ■ k(p b ■ k — ml), 

h[\ } = Amf cot ftPb ■ k(p t -k-pb-pt) - Ami tan PPt • kpb • Pt, 

h[ S 2 = -Am b m t cot f3p t ■ kp b ■ p t + Am b m t tan f3p b ■ k(p t ■ k - p b ■ p t ), 

= Ami cot [3 (2p b ■ k - p b - = p t ) - Ami t&n [3 (ml + p t ■ k), 

h!^ = — Am b m t cot /3(m 2 +p t ■ k) + Am b m t tan/3(2p& • k — p b ■ Pt), 

hf = 2m t cot(3(2p b ■ kp b ■ p t - m 2 b p t ■ k - 2(p b ■ p t ) 2 ) + 2m 2 b m t tun (3(p b ■ k — 2p b ■ p t ), 

hf = 2m 2 t m b cot (3(p b ■ k - 2p b ■ p t ) + 2m b tan(3(2p b • kp b ■ p t - m 2 b p t ■ k - 2(p b ■ p t ) 2 ), 

= 2m\ cot (3(p b ■ k — 2p b ■ p t ) + 2mlm t tan (3(p t ■ k — 2ml), 

= 2m 2 t m b cot f3(p t ■ k — 2m 2 ) + 2m b m 2 t tan/3(pf, • k — 2p b ■ p t ), 

hfp = —Am t cot (3{vr? t p b ■ k + 2p b ■ kp t ■ k) — Am\m t tan (3p t ■ k 

= —Am 2 t m b cot (3p t ■ k — Am b tan f3(mlp b ■ k + 2p b ■ kp t ■ k), 

hf = Ami cot fiPb ■ k(p b ■ k - p b ■ p t ) - Ami tan PPb • PtPt • k, 

h'd = -Am b m t cot (3p b ■ p t p t ■ k + Am b m t tan (3p b ■ k(p b ■ k - p b ■ p t ), 

hf} = Ami c °t flpb ■ k(p t ■ k — mf) — Amlm 2 t tan (3p t ■ k, 

hi2 = —Am b m\ cot (3p t ■ k + Am b m t tan (3p b ■ k(p t ■ k — ml). 
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Figure 1: Feynman diagrams contributing to the 0(a ew m^ b - ) /m'^ v ) Yukawa corrections to 
gb — > tH~: (a) and (6) are the tree level diagrams; (c) and (e) are gqq(q = b,t) vertex 
diagrams; (g) and (i) are self-energy diagrams; (k)-(m) and (o)-(q) are gbH~ vertex; (s)-(w) 
are box diagrams; (d), (/), (/i), (j), (n) and (r) are counterterm diagrams. The dashed lines 
represent H, h, A, H^, G°, for diagrams (c), (e), (g) and (i), and H, h, A, G° for diagrams 
(k), (o), (i) and (v). For diagrams (7), (m), (p), (q), (s) and (u), the dashed line (2) represents 
H and /i when the dashed line (1) is H~ , and H, h and A when the line (1) is G~ . 



t t(b) 




(d) 



Figure 2: Self-energy Feynman diagrams contributing to renormalization constants. The 
dashed line represents H,h,A,H ± ,G°,G ± in (a). 




Figure 3: The tree-level total cross sections (a) and relative one-loop Yukawa correc- 
tions (b) versus m#± at the Tevatron with y/s = 2 TeV. The solid, dashed and dotted 
lines correspond to tan/5 = 2, 10 and 30, respectively. 




Figure 4: The tree-level total cross sections (a) and relative one-loop Yukawa correc- 
tions (b) versus m#± at the LHC with y/s = 14 TeV. The solid, dashed and dotted 
lines correspond to tan/5 = 2, 10 and 30, respectively. 



